Trapped Phase-Segregated Bose-Fermi Mixtures and Their Collective Excitations 
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Recent progress in the field of ultracold gases has allowed the creation of phase- segregated Bose- 
Fermi systems. We present a theoretical study of their collective excitations at zero temperature. 
As the fraction of fermion to boson particle number increases, the collective mode frequencies take 
values between those for a fully bosonic and those for a fully fermionic cloud, with damping in 
the intermediate region. This damping is caused by fermions which are resonantly driven at the 
interface. 
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I. INTRODUCTION 

Since their first experimental observation [l[ , degener- 
ate boson-fermion (BF) systems have been realized using 
more than five different particle mixtures Q. The one 
composed of fermionic 40 K and bosonic 87 Rb particles 
is the most thoroughly studied: experiments have been 
performed on its collective excitations [|[ and on collapse 
and expansion dynamics 0, Q • Using a Feshbach reso- 
nance to tune the interspecies interaction towards repul- 
sion produced strong signs of phase segregation @, Q; 
the presence of such resonances has also been observed 
in BF mixtures of different species Despite these ad- 
vances, many new physical phenomena still await their 
experimental realization. In this context, it has been 
theoretically predicted that, for homogeneous systems, a 
multitude of phases appear on varying the interaction pa- 
rameters In trapped systems, this gives rise to states 
with a large variety of spatial structures [10|] . 

As is known from experiments with single-component 
Bose and Fermi systems, much physical information can 
be extracted by probing collective motions. Theoreti- 
cally, the collective excitations of BF mixtures have been 
intensively studied for the phase-mixed BF system [TU ] . 
inclu ding the regime close to phase-demixing [l3| . In 
Ref. [l2|, the collective excitations of phase-segregated 
Bose-Fermi mixtures across phase demixing were consid- 
ered in case of very small fermionic particle number. 

In this work, we focus on the collective motions of 
isotropically trapped, fully-segregated 40 K- 87 Rb mix- 
tures. More specifically, we investigate the eigenfre- 
quencies as a function of the fraction of the total par- 
ticles that is fermionic, a parameter which is accessi- 
ble in experiments. Using an exact analytic solution 
for the Boltzmann-Vlasov equation of a trapped, colli- 
sionless Fermi gas, we investigate monopole modes and 
their (collisionlcss) damping. To qualitatively investigate 
lower symmetry (multipole) modes we use a fully hydro- 
dynamic description. This leads to the prediction of in- 
phase and out-of-phase modes, with similar oscillations 
of inner and outer boundaries. Whereas the breathing 
modes frequencies smoothly cross over from the purely 



fermionic to the purely BEC spectrum, we find that the 
presence of a small fraction of fermions has a dominant 
effect on all other modes. 

We consider in the following Nb bosons and Np 
fermions, trapped by isotropic harmonic confinement 
Ui(r) = rriiiv?r 2 /2 with i = B, F; we shall take raflU^ = 
rriFUip such that both particle species feel the same ex- 
ternal potential, as is the case for the experiments on 
40 K and 87 Rb mixtures [1, 0, 1, [|]. Within mean-field 
theory the equilibrium state of a BF mixture with chem- 
ical potentials [Lb and jip, is fully characterized by the 
wave function ips of the Bose-Einstein condensate (BEC) 
and by the fermionic density rip. The associated grand 
potential is 0, [l(| : 
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with Gpp = ti 2 {6ir 2 ) 2 / 3 /2m Fl G B f = 2irti 2 a B F(m F 1 + 
m'g 1 ) and Gbb — ^nf? o.bb I m B the interspecies and 
bosonic intraspecies coupling constant respectively. Min- 
imization of this grand potential with respect to ni?(r) 
and iPb{?) gives the Thomas- Fermi (TF) and Gross- 
Pitaevskii (GP) equations. When, however, the parti- 
cle number is sufficiently large, the ground state config- 
uration is well described using a local-density approx- 
imation. Therefore, each point r in the trap can be 
treated locally, using the effective chemical potentials 
/Xi(r) = Hi — Ui(r), where Hi = mojfRf/2 the chemical 
potential at the center of the trap and i = B, F. 

Full BF phase segregation appears for sufficiently large 
interspecies repulsion; more specifically, it occurs when 
Gbf is larger than Gbb^f / Hb at the BF overlap re- 
gion. For all the 40 K- 87 Rb trap configurations consid- 
ered here, this requires the interspecies scattering length 
asF to be larger than 10 3 Bohr, which is well below the 
experimentally reported values of Ref. The phase- 
segregated BF trap then consists of a purely bosonic core 
surrounded by a purely fermionic shell; in between there 
is a spherical shell where the bosons and fermions over- 
lap. The thickness of this shell is generally very thin, 
of the order of the BEC healing length H/(2mB^B) 1 ^ 2 ■ 
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FIG. 1: (Left) Cross section of trap. Fermions and bosons 
are separated by the interface at radial position f. Fermions 
follow trajectories along the dashed line and reflect off the 
interface. (Right) Boson and fermion densities ub and uf 
against radial position. 



For almost all our further considered configurations this 
length is much smaller than the trapping length scales 
h/niiLOi with i — B,F. 

The influence of a thin interface on the equilibrium 
configuration can be incorporated through an interface 
tension a, an exact expression for which is given in 
Eq. (|14[) in the appendix. Mechanical equilibrium at the 
interface implies that the pressures Pi on either side of 
the interface are related by Laplace's formula: 
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Here, £ is the radial position of the interface (see Fig.[lJ. 
Finally the particle numbers of the BEC core and fermion 
shell are: 
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where K(x) = xy/l — x 2 (lAx 2 — 3 — 8x 4 ) + 3 arcsin(x). 

The GP and TF equations, together with the number 
equations ([2]) and the condition for mechanical equilib- 
rium at the interface |T]) allow us then to fully determine 
the equilibrium trap configuration. 

Having established the equilibrium configuration of the 
BF system, we now proceed to the description of its dy- 
namics. We start by discussing the dynamics of a BEC, 
which, as is well known, is governed by the equations of 
hydrodynamics. Then, both collisionless and hydrody- 
namical behavior is considered for the outer Fermi shell. 
While a collisionless treatment is experimentally more 
relevant, the succeeding hydrodynamic consideration al- 
lows to bring out some important qualitative aspects of 
the system's behavior. By combining the fermionic and 
bosonic dynamics by means of the appropriate boundary 
conditions, we then extract collective mode frequencies. 



II. DYNAMICS OF THE BEC 

As is usual, collective excitations are treated in a hy- 
drodynamical description by considering small-amplitude 
density oscillations 5ns which linearize the continuity, 



d t 5riB + V ■ {ub^b) — 0, and Euler, msn_B<9 t VB — 
VSPb — ub'S/Ub, equations. Here, tib — \iPb\ 2 is 
the BEC condensate density. Using standard meth- 
ods [3 EH* the chemical potential oscillation S/ib with 
temporal frequency u> then satisfies: 

m B nBio 2 Sfj,B + Ms V • (n B VSfi B ) = 0, (3) 

the solution to which is S/ib oc Y™r £ F(a + , a~ , I + 
3/2, [r/Rs) 2 ), with F the hypergeometric function, Y™ 
the spherical harmonic, 2a ± = z± [z 2 + 2(uj 2 /uj 2 b — i)] 1 ^ 2 , 
z = £ + 1/2 + 1 and £ a positive integer or zero. 

Here, we must impose a boundary condition at the BF 
interface. Since we are assuming that the oscillations 
occur at local thermodynamic equilibrium [2l| . we con- 
clude that the two fluids are completely immiscible and 
the interface impermeable. This immediately leads to the 
condition 



(4) 



in which is a unit radial vector normal to the interface, 
and the dot denotes the derivative with respect to time. 
That is, the interface velocity must be equal to the BEC 
velocity at the interface, ensuring vanishing of the flux of 
bosons through the interface. 



III. COLLISIONLESS FERMI GAS 

If the fermionic particles are collisionless, as expected 
at low temperatures, the appropriate dynamical descrip- 
tion is the Boltzmann-Vlasov equation with a vanishing 
collision integral: 



dtf 



; ,r ■ V v / = 0. 



(5) 



The mean- field interaction term is also absent, again 
because fermionic interactions are suppressed by Pauli 
blocking. Since collisions do not act to restore local equi- 
librium during the oscillation, the Fermi sphere is de- 
formed during the oscillation. We therefore write / = 
h+&{\p\-Pf)-v{r,x)e~ lult , where p f (r) = prnF/^r)] 1 / 2 
is the local Fermi momentum and \ — cos 6 with 9 the 
angle between r and p (see Fig.l). This parametrization 
is the most convenient for studying spherically symmetric 
oscillations. Eq. ([5]) then becomes 

— iuuv + ujfPX drV + ^f(1 — X 2 )9i r )d x L/ = 0, (6) 

with g(r) = p/r — r/p and p 2 = R 2 F — r 2 . The general 
solution to Eq. © can readily be found: 

v{r, X ) - C[r 2 (R 2 F - r 2 )(l - X 2 )V" t/2 , (7) 

with C[x] an arbitrary function of the variable x [22l |. and 

r(r,x) = arctan[2x/g(r)]/wF. (8) 

Physically, t(£,x) is the time for a particle, leaving the 
interface at an angle 9 = arccos(x) and speed p/(C)/tof, 
to return to it. 
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To fix the arbitrary function C in the general solution 
of Eq. ([7]), we must once more impose the appropriate 
boundary condition. As already explained, the interface 
is completely impermeable to both particle types; inci- 
dent fermions are therefore specularly reflected [l5|. A 
boundary condition encoding this process may be derived 
using the following arguments, after Bekarevich and Kha- 
latnikov (l6| : first, in a frame moving with the interface, 
the energy of the incident and reflected excitations must 
be equal; second, the influx of particles in this frame 
must be equal to their outflux. These two conditions 
result in 11811: 



(9) 



where ( is the velocity of the interface. One may un- 
derstand Eq. (J9j) as follows: for a stationary interface, 
C = 0, all particles incident at angle 9 will be reflected 
such that v(x) — v {~x)- If, on the other hand, the inter- 
face is moving with ( > 0, then v{x) — K~x) IS positive 
and proportional to x> that is, a) the reflected state at 
any given angle becomes more highly occupied than the 
incident state at the same angle, and b) the difference 
in occupancy is largest when % = 1, that is, for particles 
incident normal to the interface. Thus, there is an overall 
probability transfer to states closer to normal (x = 1). 
This is due to the momentum transfer from the interface 
to the reflected particles. At r = (, the solution of the 
Boltzmann-Vlasov Eq. (J7]) which satisfies condition (J9j> 



v (C, x) = m F x( [1 - i cot(W/2)] 



(10) 



Demanding local mechanical equilibrium at the inter- 
face results in the linearized Laplace equation with the 
role of the fermionic pressure fluctuation played by the 
radial component of the momentum flux tensor <5IT^, = 
/d 3 p( X p) 2 (/-/o)/[m F (27rft) 3 ]. Thus, 



5Ui 



aS((£ -!)(£ + 2) /( 2 



(11) 



- 3a5(mLO F / ' \x B + 3adfi B /(p, B () ~ 5(9 r (Pb - Pf) , 

with all quantities evaluated at the interface and 6( oc 
Y" 1 (9, <f>) the departure of the interface from its equilib- 
rium position. Solving Eqs. (fTTj) and ([4J for u> fixes the 
eigenfrequencies . 



IV. DAMPING 

Despite the absence of collisions in the fermionic gas, 
damping of the collective oscillation may occur, induced 
by the interaction of the fermions with the interface. To 
understand its origin, consider the two time scales rel- 
evant in the collisionless fermionic phase: the period of 
oscillation 2tt/uj, and t, the time of flight for a particle 
projected at angle 9 — arccos(x) from the interface (see 
Eq. ©). From Eq. HU]), v(C,x) diverges if: 

(12) 
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T i(,x) = 2im/uj, 



FIG. 2: (a) The full lines outline the region of possible values 
of r (see Eq.((S}) as a function of the interface position (. 
(b) For any multipole mode, damping will occur for values of 
5R(u;) and ( in the shaded region. The dashed lines correspond 
to the monopole modes depicted in Fig. [3] with full lines. 



with n an integer; that is, when there exist particles 
which depart from and return to the interface in exactly 
an integer number of oscillation periods. The resulting 
pole on the integration path in the pressure tensor 511^. 
and thus also in the dispersion relation Eq. (jlip , must be 
avoided by some definite prescription [17| . The appropri- 
ate manner is by deforming the contour of integration for 
STl^r so as t o pass below the pole [23| . This gives rise to 
an additional, imaginary term due to the residue in 5Il^ r ; 
in this way, the solution for co to the dispersion relation 
acquires a negative imaginary part: damping occurs. 

Can one intuitively understand when damping occurs? 
In other words, under which conditions are there particles 
which satisfy Eq. (|12p ? As we will see, whether damping 
occurs or not depends solely on the interface position and 
the real part of to, 5ft(w), and this remains true for the 
multipole modes. 

For fixed interface position (/Rf, the time r takes 
values only within a restricted interval (Fig. (2^l). Con- 
sider first radially moving particles (x = 1); for van- 
ishing Q/ Rf, they take half a trapping period to return 
to the interface, i.e., r(0, 1) = ir/tup and, as £ is in- 
creased, t smoothly vanishes as the interface approaches 
the trap edge. For (/Rf > l/\/2, particles projected 
parallel to the interface (x — 0) experience a "centrifu- 
gal force" which is insufficient to overcome the harmonic 
trapping force such that t((, 0) = 0; for general x, one 
finds < t(C,x) < t(C, 1). Conversely, if (/R F < 1/V2, 
a particle with x = performs exactly half an orbit 
between two collisions with the interface and therefore 
t(£,0) = tt/ujf', in this case, t(£, 1) < t((,x) < tt/ujf 
for general x- 

From the possible values of r, it is straightforward to 
obtain the values of uj and £ for which damping occurs 
(see shaded region in Fig.J^b). For small (/Rf, we have 
seen that r » it/ujf; but this is exactly half the period of 
a collisionless Fermi gas. Damping is therefore expected 
for the lowest monopole mode when N B <C Np- On 
the other hand, when (/Rf ~ 1, due to a vanishing 
r, Eq. (TT2)) is satisfied only for large values of u. This 
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coiiisioniess Fermi Fraction Q=N F /(N B +N F ) 

Damped Coiiisioniess — — • Hydrodynamic 

FIG. 3: The frequencies u of the 1 = monopole modes as 
a function of the Fermi fraction Q = Nf/{Nb + Nf). The 
full black and grey lines correspond to damped and undamped 
monopole modes, respectively, both calculated using collision- 
less fermion dynamics. The inset shows the damping rate T 
of the second monopole mode. Dashed lines are calculated 
using a hydrodynamic Fermi shell. 
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FIG. 4: The frequencies u> of the multipole (£ 7^ 0) modes as 
a function of the Fermi fraction Q = Nf/(Nb + Nf)- All are 
found using a fully hydrodynamical trap and involve in-phase 
(full lines) and out-of-phase (dashed lines) motions of inner 
and outer boundary. 



VI. OUT-OF-PHASE MODES 



inhibits damping for the lowest modes in case Np <C Nb ■ 
The damping we observe is somewhat analogous to 
Landau damping, which corresponds to energy trans- 
fer from collective modes into single-particle incoherent 
modes. In our case, no clear distinction exists between 
collective and single-particle modes in the Fermi gas and 
energy goes into incoherent excitations in the coiiisioniess 
gas. 



V. RESULTS 

We consider trap configurations containing 10 6 parti- 
cles and characterized by the parameters asB = 5.23 
nm and ujb = 570 Hz [2J|. In Fig. [3J we show with 
full lines the real part of the frequencies of the two low- 
est monopole modes as a function of the Fermi fraction 
Q = Nf/{Nb + Np). The grey and black lines indicate 
the damped and undamped frequencies respectively. The 
lowest monopole mode smoothly crosses over between the 
value y/Eu>B for a fully superfluid bosonic trap at Q = 0, 
to the value 2u>p of the fully coiiisioniess (fermionic) trap 
at Q = 1 [l4j]- As explained earlier, the lowest monopole 
mode is damped when C/Rf » 0; indeed, for Q > 0.85, 
very weak damping is present with a maximum damping 
rate of about 10 _3 cj_b. As expected, the second monopole 
mode frequencygoes over from ^/\Aujb at Q — to 
Aujf at Q = 1 14]. Its variation, however, is nonmono- 
tonic and considerable damping occurs in the interval 
0.02 < Q < 0.4; the associated damping rate T is shown 
in the inset of Fig. [3J The absence of damping for low val- 
ues of Q is in full agreement with the explanation given 
above. In Fig. [2b, we plot the evolution of u> as a function 
of £ (dashed lines) for both monopole modes; damping 
occurs in the shaded region. 



While a coiiisioniess treatment is experimentally more 
relevant, the following hydrodynamic considerations for 
the fermions allow us to bring out some important qual- 
itative aspects of the system's dynamical behavior. 

The bulk solutions to the continuity and Euler equa- 
tions for fermions are similar to those for the BEC [la |. 
The Laplace Eq. used for the coiiisioniess case still 
holds but with <5IT^, replaced by ripS/iF and the imper- 
meability of the interface implies <5£ = • Vp^. 

In Fig. [3j we show with dashed lines the frequencies 
for the breathing modes t = using a hydrodynamical 
Fermi shell. In Fig. [H on the other hand, the lowest 
multipole (£ ^ 0) modes are shown. Both the multipole 
modes in which the interface and outer boundary move 
in phase (full lines) and out of phase (dashed lines), are 
present [T^|. Both the IP and OOP modes have in-phase 
velocities at the BF interface and are both generic to a 
two-component system [18l . Il9j . An interesting result is 
that the frequencies of all the £ ^ modes vary strongly 
at small Q and remain approximately constant at larger 
Q; this in contrast with the monopole modes. 



VII. DISCUSSION 

The mass difference between the 40 K and 87 Rb par- 
ticles may cause a gravitational sag to be present. In- 
corporating such effect is beyond the scope of this work. 
However, we expect gravity not to affect the validity of 
our qualitative results such as the damping mechanism 
and the presence of out-of-phase modes. Note also that 
the gravitational sag may be quenched experimentally 
by counterbalancing external potentials or by tight ra- 
dial confinement. 

A finite temperature will change the collective mode 
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frequencies mostly through a change in the fermion dy- 
namics. We expect 5R(w) to be changed little with 
temperature, while damping is expected to occur for a 
broader range of Q and occurs by the same mechanism 
as described here [201 ] 

Our approach fails when the interface is as thin as the 
outer fermion shell; we find that this happens for Fermi 
fractions below 2 x 10~ 3 . In that case, damping of the 
collective excitations happens by usual Landau damping 
throughout the BF overlap region [l2j] . Our work is com- 
plementary to Ref. [l2j in the sense that they study the 
process of phase demixing for a Fermi fraction 4 x 10~ 5 
whereas we focus on the variation of the Fermi fraction 
in case of complete phase segregation. 

VIII. CONCLUSION 

We present monopole mode frequencies for the trapped 
phase-segregated 40 K- 87 Rb mixture. A damping mecha- 
nism is identified which is generic to all collective modes. 
Within a fully hydrodynamic approximation we then ob- 
tain mode frequencies for modes involving out-of-phase 
and in-phase motions of inner and outer boundaries. All 
our results are only very weakly dependent on the surface 
tension. 
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X. APPENDIX 

In order to obtain the surface tension of a BF in- 
terface, one may consider an infinite space containing 



fermions at x — > — oo and bosons at x — > +00, bound 
by a flat interface which is parallel to the y — z plane. 
Both species are taken to be untrapped (Ub = Uf = 
and at fixed chemical potentials \xb and \xf- The ap- 
propriate boundary conditions for tpB(x) and np(x) are 
V>b(-oo) = n F (+00) = while ip B (+°o) = \Z^b/G B b 
and np(— 00) = (^f/Gff) 3 ^ 2 - We define the BEC co- 
herence length £b = ^/[^.ttib^b] 1 ^ 2 an d the dimension- 
less parameters k 2 = Gbb^f/Gbf^b and (/)b(x) = 
iPb(x)/ ' \J ^b/Gbb- Using these variables, is useful to 
introduce the first integral associated with the GP and 
TF equations for our system under consideration: 

2e B (d<j> B /dx) 2 = (1 - 4>lf - G(k, </> B ), (13) 

where £(k,<M = (1 - <p 2 B / k 2 ) 5 / 2 0(k - (j) B ) and 9 the 
Heaviside function. Using Eq. (fl"3|) , the interface tension 
which is the excess grand potential per unit area of the 
interface can be recast in the form: 



a=4P£! / dxidfa/dx) 2 , 



Using the GP and TF equations and the transformation 
dx = {d4>B / dx)~ l &4>B, the interface tension a associated 
with a flat BF interface is: 

ct(k) = 2V2P£ B J d^a-^!) 2 -^,^), (14) 

The surface tension is minimal at K = 1 (a = 0.52P£ B ) 
and reaches its maximal value 4v / 2P£s /3 when k = 0. 
Note that one may use this exact result for the surface 
tension in a trap when the chemical potentials [Ib and 
fip are evaluated at the interface. 
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